Electronic and transport properties of rectangular graphene macromolecules and 

zigzag carbon nanotubes of finite length 
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We study one dimensional (ID) carbon ribbons with the armchair edges and the zigzag carbon 
nanotubes and their counterparts with finite length (OD) in the framework of the Hiickel model. We 
prove that a ID carbon ribbon is metallic if its width (the number of carbon rings) is equal to 2 + 3n. 
We show that the dispersion law (electron band energy) of a ID metallic ribbon or a ID metallic 
carbon nanotube has a universal sin-like dependence at the Fermi energy which is independent of 
its width. We find that in case of metallic graphene ribbons of finite length (rectangular graphene 
macromolecules) or nanotubes of finite length the discrete energy spectrum in the vicinity of e = 
(Fermi energy) can be obtained exactly by selecting levels from the same dispersion law. In case 
of a semiconducting graphene macromolecule or a semiconducting nanotube of finite length the 
positions of energy levels around the energy gap can be approximated with a good accuracy. The 
electron spectrum of OD carbon structures often include additional states at energy e = 0, which 
are localized on zigzag edges and do not contribute to the volume conductivity. 

PACS numbers: 73.22.-f, 73.20.-r, 73.23.-b 



Carbon based materials of nano-size in the form of 
carbon nanotubes (CNTs) are known for several years 
and attracted much attention of researchers because of 
their unusual electronic properties [3, [3, [I]- Recently, 
progress in the fabrication of other graphene-based lower 
dimensional structures has been reported 0, [f|. This 
put forward such nano-scaled quantum objects as one- 
dimensional (ID) carbon ribbons (CR) MM !, H 03 

SIEf. 




and zero-dimension (OD) carbon dots 
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Graphene - a two-dimensional (2D) carbon material - 
was first isolated by micromechanical cleavage of graphite 
Its planar hexagonal lattice is formed by sp 2 hy- 
bridized carbon bonds. Although graphene is the build- 
ing block of many carbon allotropes, its electronic struc- 
ture differs from other carbon materials. At present, the 
electronic structure and transport properties of ID CNTs 
are well understood theoretically [13| and the focus is 
shifted towards CRs and OD carbon objects In par- 
ticular, it is known that in the framework of the tight- 
binding model both CRs and CNTs (ID nano-materials), 
can be either semiconducting with a size dependent gap 
or metallic It should be also noted that recent 

density functional calculations within local density ap- 
proximation (DFT-LDA) predict that all armchair CRs 
are semiconducting, with one group showing small en- 
ergy gaps [H, [l?], EH ■ We will not discuss this issue here 
and limit our consideration by the tight-binding (ID) 
and Hiickel model (OD). There, the rule of metallicity 
for armchair CRs was formulated under assumption that 
the ribbons are wide enough so that one can use solu- 
tions obtained for graphene semiplane [l5|, El . Below we 
refine the procedure and obtain the electronic solutions 
which are equally applicable to CRs of small and large 
width. Furthermore, we show that for the chosen class 
of ID metallic carbon systems the dispersion law of the 
electronic band crossing the Fermi level can be obtained 




FIG. 1: Boundary conditions for a ID carbon ribbon (A) and 
the same ribbon placed in the graphene sheet (B). Carbon 
sites on lines L\ and L2 have zero coefficients of expansion. 



analytically [see Eq. © below]. Later we generalize the 
law for OD systems [Eq. (|14a|) . (|14b|) below]. Throughout 
the letter we employ the Hiickel model and limit our- 
selves to the class of zigzag CNTs and CRs with the 
armchair profile. These two materials are closely related 
with each other. Indeed, by rolling up an ID armchair 
CR one obtains a ID zigzag CNT. (Terms armchair for 
the ribbon and zigzag for the nanotube are confusing here 
since they apply to different characteristics: armchair - 
to the edges of the ribbon, while zigzag - to the circum- 
ference of the nanotube.) On the other hand, all CRs of 
finite length can be considered as rectangular graphene 
macromolecules (RGM) whose electronic properties are 
important for designing various nano-materials. 

The energy spectrum of nanographite materials (ID 
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CRs, ID CNTs, OD RGM, OD CNTs) can be obtained 
from the dispersion relation of graphene [201 ] : 



SF 



± /3| e -*<»*»V^/2 +2008 0^/21, 



(1) 



where a = \/3dcCt and dec is the carbon-carbon dis- 
tance. Here the 2ir factor is incorporated in k, and —/3 
stands for the Hiickel transfer integral (or — t in the tight- 
binding model), so that (3 > 0. At the K~ point of the 
Brillouin zone of graphene two bands intersect: 
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(2a) 
(2b) 



We start with studying energy spectrum of ID CRs. 
The problem is how to reduce the problem to that for 
graphene. In the unit cell of a zigzag carbon ribbon 
there are four carbon atoms with 2 nearest neighbors 
while in graphene each carbon atom has 3 nearest neigh- 
bors, Fig. [TJ Therefore, if we consider the ribbon as a 
part of the graphene 2D plane, Fig.Q}3, the equations for 
the edge atoms of the ribbon are modified. For example, 
for site in Fig. [TJ3 we have: 



-c = -(ci + c 2 + c 3 ), 



(3) 



where Cj, i = 0, 1, 2, 3 are coefficients of the expansion of 
graphene wave function ip g . However, if we put C3 = 0, 
then Eq. will be identical to that for the ribbon. Thus, 
we can select the solutions for the CR out of the graphene 
solutions by requiting cj — for all carbon sites j on lines 
L\ and Li and by repeating the resulting pattern in the 
x— direction. The lines with tp r — completely separates 
neighboring nanoribbons from each other, since there is 
no interaction between them. 

We arrive at the following conditions of quantization: 



cos{Dk' x ) = 0, 
sm{Dk' x ) = 0. 



(4a) 
(4b) 



Here D = md (m is an integer) and d = y/3dcc/2- (Both 
distances are shown in Fig. [TJ3.) The number of carbon 
rings in the x— direction is M = m — 1. It is clear that 
the CR will be metallic if line (k' x , k' y ) in k— space goes 
through the if-point, Eqs. (J^aJ), pfe]). Solving Eqs. (|Iaj) . 
(I4bl) leads to 



I lb 

ak„ = 2n — , 
m 



(5) 



where n is an integer number. Condition ([5]) coincides 
with Eq. (gal), (E9 only if 



(6a) 
(6b) 



m = 3m' , 
71 = 2m'. 

This immediately gives 

Af = 3m'-1. 



(7) 




FIG. 2: Unit cell of ID metallic carbon ribbon. Plus signs 
(sites 6, 10) and minus signs (2, 8) refer to eigenvector Vi at 

EF = 0. 




FIG. 3: Unit cell of ID metallic carbon ribbon containing 
8 carbon rings. Plus and minus signs schematically indicate 
one of two eigenvectors at ef = 0. 



Eq. ([7]) determines which CR is metallic. We get: J\f=2, 
5, 8,... . Substituting (k' x = fcf ,k' y = k) in Eq. ^} leads 
to the dispersion relations for two bands which cross at 
the Fermi energy: 



= Ef ± 2(1 sin ^j— , 



(8) 



where X = 3dcc (X is the modulus of the basis vector 
in the y— direction). The corresponding density of states 
(DOS) per unit cell (for both spins projections) in the 
neighborhood of sf = is given by 
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(9) 



Eq. © coincides with that for metallic CNTs [13| . 

It is instructive to study the two electronic states at 
the Fermi energy, when kp = 0. We start by considering 
the simplest possible case: N — 2, Fig. [2l Then the 
following explicit eigenvectors can be found: 

Vx = \ {0,-1,0,0,0,1,0,-1,0,1}, (10a) 

V 2 = ~ {-1,0,1,0,-1,0,1,0,0,0}. (10b) 

V\ is visualized in Fig. [T] by putting plus and minus signs 
at corresponding carbon sites. V2 is obtained from V\ 
through the mirror reflection at the xz plane, which fol- 
lows from the symmetry of the electronic system. In this 
way one can construct two eigenvalues for any metallic 
CR. As an example in Fig. [3] we schematically draw one 
eigenvector for a CR with the 8 ring unit cell. Notice, 
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that all three blocks have the structure of V\ , Fig. [2] The 
blocks are connected via chains of carbon sites - lines L\ 
and 1/2 - where Cj — 0. As before, the second eigenvector 
is obtained from the first through the xz— mirror reflec- 
tion. Thus, we can build two eigenvectors only if number 
of blocks is Af = 2 + 3n, which also gives the rule ((7|). 

From Eq. ([1]) one can derive the expression for energy 
gap EgiiM^ 3m' — 1. There are two cases:(l) Af = 3m' 
and m = 3m' + 1, and (2) Af = 3m' — 2, m = 3m' — 1. 
In both cases the energy gap is given by 



E g = 20\l 



cos — Vosm — | 



(11a) 



The value of 5<j) depends on which branch - Eq. (|4a[) or 
Eq. (|4b[) - is considered. [It is a measure of deviation from 
the if— point, Eq. (J^ajb) in the zone-folding model.] For 
cos-like branches, Eq. [[la]), we get n = 2m', Eq. (l6b|) . 
and 

(lib) 



rtos '3(jV+l)' 

Here the first choice of Af corresponds to minus sign, 

while the second choice of Af to plus sign. For sm-like 
branches, Eq. (|4aj) . we get 



*>sin = 40 co;s , n = 2m' 
*> sin = 2<f) cos , n = 2m' ± 1 



(11c) 
(lid) 



By comparing Eq. (|llb|) with Eq. (fTTc|) and Eq. (|lld|) 

we conclude that the gap is due to cos-like branches, 
Eq. (jllbp . If S(f) cos <C 1, which is often the case, 



E g = (3 



4tt 1 

VsAf + T' 



(12) 



Thus, E g ~l/AT as it was the case for CNTs [ijf. 

The rule for the metallicity of armchair CRs, Eq. Q, 
is very different from the the rule for the corresponding 
zigzag CNTs characterized by the pair of indices [Af, 0] : 



Af = 3m'. 



(13) 



It is clear that the rules do not overlap meaning that if 
one rolls up a nanotube from a metallic ribbon, then 
the resultant nanotube will not be metallic and vice 
versa. This conclusion deserves a more detailed explana- 
tion. First, we notice that the metallicity rule for CNTs, 
Eq. (fl~3"j) , is obtained from the cyclic condition in the 
a;— direction which allows for sin- and cos- like depen- 
dencies, while in case of ribbon only sin- like functions 
are allowed. [One can prove that Eq. (|4a[) does not in- 
clude the if— point, Eq. (j2aj) .] The sin-dependence in 
the x— direction implies opposite signs for coefficients Cj 
belonging to edge sites j on opposite sides of CR. By 
rolling up a nanotube the carbon sites on opposite edges 
of ribbon should coincide, which in turn destroys the odd 
solutions. The even solutions satisfying the cyclic condi- 
tion survive the rolling up procedure but none of them 
has energy at ep — 0. 



TABLE I: Length dependence of the energy spectrum of 
metallic RGMs (Af = 2, 5, ...) and CNTs (Af = 3, 6, ...). 



HOMO- 



length, in dec 
299 149 29 



1 


-0.00783 


-0.01563 


-0.07661 


2 


-0.02350 


-0.04689 


-0.22937 


3 


-0.03917 


-0.07813 


-0.38078 


4 


-0.05483 


-0.10935 


-0.52996 


5 


-0.07049 


-0.14055 


-0.67603 


6 


-0.08615 


-0.17172 


-0.81814 


7 


-0.10180 


-0.20284 


-0.95544 



We now turn to 0D objects - carbon macromolecules 
and nanotubes of finite length. Unlike ID electron sys- 
tems characterized by electron energy band structure, 
they have discrete energy spectra. However, these spec- 
tra are closely related with electron bands which we have 
already considered. In particular, energy levels of RGMs 
and CNTs near the Fermi energy are described by the 
following expression: 



where 



i + _(n) = ±2/3 sin 



k(n) — 8k In + 



Xk(n) 



(14a) 



(14b) 



n = 0, 1, 2,... and 5k = ir/L, and L = Lq + Idccl^- Here 
Lq is the length of nanotube or ribbon (maximal distance 
between carbon atoms in the y— direction, which has the 
armchair profile). We want to stress that Eqs. (|14a[) . 
(I14bp are exact, see Appendix. The electron spectrum 
given by (|14a|) is independent of width which is consis- 
tent with the situation observed for ID objects. It is also 
worth noting that the spectrum of 0D carbon objects con- 
sists of many discrete levels and instead of Fermi level we 
should speak of highest occupied (HOMO) and lowest un- 
occupied (LUMO) molecular orbitals. This implies that 
formally we can not speak of metallicity and semiconduc- 
tivity. However, in Tables Hl lIIII we retain these terms in 
a loose sense, because from one side, it shows relations 
with corresponding ID objects and from the other side, 
the existence or nonexistence of a large energy gap at 
the HOMO-LUMO region remains one of the important 
characteristics of these systems. 

It is clear that Eqs. (|14ap and (I14b[) can be considered 
as a discretization of ([5]). Following this route we can 
derive an approximate expression for nonmetallic RGMs 
and CNTs. First, we recall that for ID systems there are 
various bands, which we have discussed already while 
calculating E g . In general their energy is 



£+.- 



3 



±\e~ 



cos ■ 



V3i 



(15) 



where 5<j> = 5<p S in or & 
for RGMs is given by 



cos 



The highest occupied band 
Eq. (jllbp . Starting with 
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TABLE II: Length dependence of the energy spectrum of 
semiconducting RGM, TV = 10. (L = 3.5d C c, 5 = 0.9). 
"Exact" refers to the strait-forward Hiickel calculations, "ap- 
prox" to the values of Eq. (llGal b) . 





299 dec 


214 dec 


HOMO- 


-i exact 


approx. 


exact 


approx. 


1 


-0.08271 


-0.08258 


-0.08371 


-0.08356 


2 


-0.08689 


-0.08691 


-0.09069 


-0.09100 


3 


-0.09349 


-0.09388 


-0.10143 


-0.10250 


4 


-0.10210 


-0.10294 


-0.11501 


-0.11687 


5 


-0.11230 


-0.11360 


-0.13064 


-0.13318 


6 


-0.12374 


-0.12544 


-0.14771 


-0.15079 


7 


-0.13611 


-0.13816 


-0.16580 


-0.16928 



TABLE III: Length dependence of the energy spectrum of 
semiconducting CNT [10,0] (Af = 10, L = 3.75d C c, 8 = 
0.9). "Exact" refers to the strait-forward Hiickel calculations, 
"approx" to the values of Eq. (|16al b). 



HOMO- 



299 dec 
exact approx. 



214 
exact 



dec 
approx. 



-0.17634 
-0.17864 
-0.18241 
-0.18756 
-0.19400 
-0.20161 
-0.21027 



■0.17623 
0.17847 
0.18226 
0.18750 
0.19406 
0.20182 
0.21065 



-0.17691 
-0.18088 
-0.18733 
-0.19602 
-0.20671 
-0.21910 
-0.23296 



-0.17673 
-0.18067 
■0.18723 
-0.19615 
-0.20712 
-0.21982 
-0.23398 



Eq. (I15|) and taking inspiration from the relation be- 
tween ((8]) and (|14ap we obtain the discretized version 
for e + ,-(n) of semiconducting RGMs or CNTs of finite 
length: 



±e 



£+,-(") 
k'{n) = Sk (n + S). 



<"■' r ^4> i 
cos — — vosm— I, (16s 



(16b) 



Here n = 0, 1, 2,... is integer, Sk = n/L, and L = L + SL. 
In fact, in Eq. (]16b[) SL ~ dec an d 6^1 are phenomeno- 
logical parameters which should be found by fitting a 
data set of e(n). The set can be taken from a Hiickel 
calculation of RGMs and CNTs or even from a more ad- 
vanced calculation (like density functional, for example). 
In Tables [H] and IIIII we compare energy spectra given by 
Eq. (|16al b) [approx.] with straight-forward Hiickel cal- 
culations [exact]. 5(f) in Eq. (|16a[) is 4> cos , Eq. (jllb|) 
graphene molecules, and 



for 



= ± 



2tt 



(17) 



in case of semiconductiong nanotubes with width TV = 
3m' ± 1. The accuracy of Eqs. (|16ap , (|16b[) is of the 
order of 10~ 3 /3 for 10 top occupied molecular orbitals (or 
HOMO-i, i = 1 - 10), Tables HI HIl 

Finally, we remark on the doubly degenerate HOMO 
level at e « [1, [TH, The nature of the states 

which are always present in the electronic spectrum of 
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FIG. 4: Conductance versus RGM length at ef = 0. 



carbon ribbons and nanotubes of finite length, is an ob- 
ject of intense research @, [H [H [HI [H, [H, ^ . These 
states are edge ones, because as follows from Eqs. (U4a[) . 
(fl4b)l ehomo-i = -2[3sm{X5k/9,) ^ 0. The coeffi- 
cients of the wave function expansion, Cj, in this case 
quickly fall to zero as we move away from the edges, and 
the states do not contribute to the conductivity in the 
y— direction, Fig. [T] To demonstrate it we have calcu- 
lated the conductance of RGMs using the Landauer for- 
mula 



25j . For a qualitative treatment, we define the self 
energy within the broadband approximation |26|, con- 
sidering it as an energy- independent imaginary constant, 
S = iA (A = 1 cV). In Fig. 0] we plot conductance as a 
function of length at sp = for RGMs of various width. 
Our results show that metallic RGMs (n = 2, 5, 8) have a 
weak dependence on length. Their conductances coincide 
starting with a rather short length of 2-6 rings. For other 
non-metallic RGMs the calculations demonstrate an ex- 
ponential decrease of conductance with length. It is also 
worth noting that the conductance of metallic RGMs is 
always the same and equal to the conductance of the 
metallic RGM with the minimal width of 2 rings. (A de- 
tailed qualitative and quantitative analysis will be given 
elsewhere.) 

In conclusion, wc have studied the electronic spectrum 
of the ID armchair CRs and zigzag CNTs and their 0D 
counterparts of finite length in the Hiickel model. We 
have found the solutions by reducing the problem to that 
for graphene with appropriate selection rules imposed 
by boundary conditions. In the vicinity of the HOMO- 
LUMO energy region (e ~ 0) we have found the exact 
expression for energy spectrum of metallic nanosystems, 
Eq. (|14al b) and approximate energy spectrum in case of 
semiconducting materials, Eq. (|16al b). Finally, we have 
calculated the conductance of some RGMs and investi- 
gated the role of edge states. 

The authors would like to thank D. S. Kosov for helpful 
discussions. 
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' \----»H+l 



FIG. 5: Effective ID chain of N carbon atoms [i = 1, 2, N] 
within the /-ribbon. Also shown are two auxiliary carbon 
atoms [i = and i = (N + 1)], see text for details. 



APPENDIX 



Here we discuss the derivation of Eq. (!14ai bV For 
a ID armchair CR or ID zigzag CNT each of the dis- 
crete values of k' x determines an electronic band. The 
metallic band is given by Eq. (f2"a|) and obtained for the 
(|4b[) condition of quantization. This condition implies 
the sin[(27r/3)(a;/c?)] modulation in the x— direction with 
zero coefficients at x/d = 3/ where I — 0,1,2,.... (Dis- 
tances and axes are shown in Fig.[T]) Each index / defines 
a carbon ribbon in the y— direction with nonzero coeffi- 
cients at x\/d = 1 + 3/ and X2/d = 2 + 3/. The lines 
with zero coefficients imply that each ribbon can be con- 
sidered as independent. In fact, they are equivalent due 
to the modulation condition, Eq. jib)) . We can use this 
property and work only with one ribbon shown in Fig. [5] 
and later reconstruct the solution for the whole system. 

Now we consider one ribbon in the y— direction which 
is equivalent to a one dimensional (ID) chain of carbon 
atoms (i — 1,2, ...,N). Considering a solution with the 
coefficients Cj of the wave function expansion, we obtain 
for them the following relations 



-(Cj-i +c j+ i) 



(A-1) 



where j = 2,3, (N — 1). The problem arises due to 
the two boundary atoms: j = 1 and j = N. To solve 
this task we will use the trick which we have already 
applied to ID ribbon with Eq. ([3]). That is, we introduce 
two auxiliary carbon atoms jbx — and j2& = (N + 1) 
and consider an infinite ID carbon chain beyond them. 
Notice, that for the following the real shape of the ID 
chain is immaterial, it can be equally thought of as a ID 
linear chain of carbon atoms. For an infinite chain the 
general solution is Cj ~ exp(—(j)' j), where <fi' = dec k' y 
and k' y is an effective wave number. Then the coefficients 
at jb should be zero, i.e. Co = cjv+i = 0. This gives a 
condition of quantization, 



sin(L y k' y ) = 0, 



(A-2) 



where L y = (N + 1) d cc . From (|A-2[) we get 



k' = — = - 
y L y r (N + l)d CC r ' 



(A-3) 



Here positive integer r — 1,2, ...,7V. The energy of the 
ID chain is obtained from Eq. (| A- 1 [1 : 



e/p = -2 cos(0') = -2 sin Q 



The latter relation can be written as 
e/P = -2 sin 



7T 1 

(n + - 

N+V 2 



where 



2n = N — 2r. 



(A-4) 



(A-5) 



(A- 



Notice that n = N/2 - 1, (N/2 - 2), -N/2. (N is 
even, because N = AN^ex, where N^ex is the number of 
hexagons in the y— direction.) First N/2 electron states 
are occupied by N 7r-electrons. Taking into account that 
the length of CR or CNT in the y— direction, 



3 -N 

4 



1 d 



cc, 



(A-7) 



Eq. (|A-5|) can be rewritten as 



e/P = +2 sin 



/ 3d C c 



V 4 (Lo + ldcc/±) 



(n' + -) ,(A 



where the solution with the minus sign refers to the oc- 
cupied states at the Fermi level: HOMO (n' — 0), and 
the HOMO-n' levels (ri = 1, 2, N/2.) The solution 
with the plus sign refers to the unoccupied states: LUMO 
(n' = 0), and the LUMO+n' levels (n' = 1, 2, N/2.). 
Eq. (|A-8|) is equivalent to Eq. (|14al bV Finally, we would 
like to notice that the equation implies a certain size re- 
lation of RGM, L y = Lq > L x . The latter condition is 
needed to assure that a set of levels around the Fermi 
level is associated with a ID metal band and separated 
from other k' x — levels (other ID bands). If it is not so 
(L y = Lo < L x ) then the levels described by (|A-8|) do 
exist but they are not grouped together. They are mixed 
up with other k' x — levels. 
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